UCB-PTH-09/18 



Cosmic Signals from the Hidden Sector 

Jeremy Mardon, Yasunori Nomura, and Jesse Thaler 

Center for Theoretical Physics, Department of Physics, 
University of California, Berkeley, CA 94720 
and 

Theoretical Physics Group, Lawrence Berkeley National Laboratory, Berkeley, CA 94720 

Abstract 

Cosmologically long-lived, composite states arise as natural dark matter candidates in 
theories with a strongly interacting hidden sector at a scale of 10 - 100 TeV. Light axion- 
like states, with masses in the 1 MeV - 10 GeV range, are also generic, and can decay 
via Higgs couplings to light standard model particles. Such a scenario is well motivated 
in the context of very low energy supersymmetry breaking, where ubiquitous cosmological 
problems associated with the gravitino are avoided. We investigate the astrophysical and 
collider signatures of this scenario, assuming that dark matter decays into the axion-like 
states via dimension six operators, and we present an illustrative model exhibiting these 
features. We conclude that the recent data from PAMELA, FERMI, and H.E.S.S. points to 
this setup as a compelling paradigm for dark matter. This has important implications for 
future diffuse gamma ray measurements and collider searches. 
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1 Introduction and Summary 

Weak scale supersymmetry is a very attractive candidate for physics beyond the standard model. 
It stabilizes the weak scale against potentially large radiative corrections, leads to successful 
gauge coupling unification, and predicts a plethora of new particles accessible at the LHC. This 
framework, however, also suffers from several generic cosmological problems, associated with 
overproduction of gravitinos or late decay of the field responsible for supersymmetry breaking [HO 
|3l H]. In many supersymmetry breaking scenarios, this requires a rather low reheating temperature 
after inflation, making it difficult to explain the observed baryon asymmetry of the universe. 

A simple way to avoid these cosmological problems is to assume that the gravitino is very 
light. If the gravitino mass satisfies 

m 3/2 <O(10eV), (1) 
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then the gravitino is in thermal equilibrium with the standard model down to the weak scale. 
The resulting gravitino abundance is small [5] and consistent with structure formation [6]. This 
allows for an arbitrarily high reheating temperature, and thus baryogenesis at high energies such 
as thermal leptogenesis [TJ. 

It is remarkable that this simple cosmological picture is precisely the one suggested by arguably 
the simplest scheme for supersymmetry breaking. Suppose that supersymmetry is dynamically 
broken at a scale 

\^F^ O(10 - 100 TeV), (2) 

giving 7713/2 = F/v^Mp! w 0(0.1 - 10 eV), where M Pi ~ 2.4 x 10 18 GeV is the reduced Planck 
scale. Then, if the sector breaking supersymmetry is charged under the standard model gauge 
group, gaugino and scalar masses of order (g 2 /167t 2 )a/F m O(100 GeV - 1 TeV) can be generated 
through standard model gauge loops [SJ |9| . The generated squark and slepton masses are flavor 
universal, and thus solve the supersymmetric flavor problem. The supersymmetric Higgs mass 
(fi term) can also be generated through direct interactions between the Higgs fields and the 
supersymmetry breaking sector [10l[TT] or by the vacuum expectation value of a singlet field [12]. 

What are the experimental signatures of this simple supersymmetry breaking scenario beyond 
its indirect implications on the superparticle spectrum? In this paper we advocate that the 
scenario may lead to distinct astrophysical and collider signatures, due to the following features 
of the dynamical supersymmetry breaking sector that may appear under rather generic conditions: 

Quasi-stable states: One immediate consequence of the present framework is that dark matter 
cannot be the lightest supersymmetric particle, which is the very light gravitino. Dark 
matter, however, can arise naturally as a (quasi-) stable state in the supersymmetry breaking 
sector [13]. Let 777dm be the mass of this state. The annihilation cross section is then naturally 
(av) ~ (l/87r)(ft 4 /m DM ), where k represents typical couplings between states in the strong 
sector. For 777dm ~ O(10 - 100 TeV), natural values for the coupling k ph 0(3 - 10) give 
(av) ~ (l/87r)(l/TeV 2 ), which leads to the correct thermal abundance for dark matter, 
Odm ^ 0.2. 

In general, a strongly interacting sector of a quantum field theory often possesses enhanced 
global symmetries, such as baryon number and flavor symmetries. These symmetries can 
lead to "stable" states if they are not broken by the strong dynamics. It is, however, quite 
possible that these symmetries are not respected by physics at (much) higher energies, such 
as at the gravitational scale. The "stable" states then become quasi-stable states, decaying 
through higher dimension operators on cosmological timescales. 

Light axion-like states: A light state appears as a pseudo Nambu-Goldstone boson when 
an approximate global symmetry is spontaneously broken in the strongly interacting sector. 
Since our strong sector is supposed to break supersymmetry, a generic argument of Ref. [H] 
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suggests that it may possess an accidental U(l) R symmetry which is dynamically broken 
(even if supersymmetry is broken in a local minimum). This then leads to a light R axion, 
whose mass is typically of 0(1 - 100 MeV) if the mass dominantly arises from a constant 
term in the superpotential canceling the vacuum energy of order F 2 [15J. 
More generally, it is not hard to imagine that the sector possesses enhanced approximate 
global symmetries. If these symmetries are spontaneously broken, light axion-like states with 
masses much smaller than the dynamical scale will appear. The masses of these states are 
then controlled by the size of explicit breaking of the corresponding symmetries. 

Couplings to the Higgs fields: If the supersymmetry breaking sector yields an R axion, it 
generically couples to, and thus mixes with, the Higgs fields. This is, in fact, the case even 
if the Higgs fields are not directly coupled to the supersymmetry breaking sector because 
the holomorphic Higgs mass-squared (Bfi term) obtains loop contributions from gaugino 
masses, which are necessarily i?-violating. For other axion-like states, mixing with the 
Higgs fields can arise if the Higgs fields are directly coupled to the supersymmetry breaking 
sector, making the corresponding symmetries Peccei-Quinn (PQ) symmetries. This is well 
motivated, since such couplings are often needed to generate the \i term. 

It is interesting to note that all the ingredients above appear in QCD, a known strongly coupled 
system in nature. For the first two, one should simply think of protons, neutrons, pions, and kaons 
(in the appropriate limits where quark masses are small or weak interactions are superweak). Even 
direct Higgs couplings do exist, although the Higgs boson is much heavier than the dynamical 
scale of QCD, while it is much lighter than the dynamical scale considered here. 

The structures described above could be manifested in various experiments. In particular, 
couplings between the Higgs and axion-like states provide a potential window to probe the su- 
persymmetry breaking sector directly. Since the sector may contain dark matter as well as light 
states, possible signatures may appear both in astrophysical and collider physics data, and in this 
paper we consider the following two classes of signatures: 

Cosmic ray signals from decaying dark matter: If dark matter is a quasi-stable state 
in the supersymmetry breaking sector, its decay may lead to various astrophysical signa- 
tures. Assuming the decay occurs through the light states, the final states can be mainly 
leptons, explaining the excess of the positron to electron ratio observed in the PAMELA 
experiment [16J, along the lines of [T71 [181 EH]- The required lifetime of order 10 26 sec is 
obtained if the decay is caused by a dimension six operator suppressed by the unification or 
gravitational scale [20J. 

There are characteristic features for this explanation of the PAMELA excess which are 
being tested in current observations. First, since the decay occurs through light states, it 
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Symmetry 


Spontaneous 
Breaking? 


Explicit 
Breaking? 


Consequences 


R or PQ 


Yes 


Sunereravitv or 
d = 5 at Mi 


Axion-like state with: 
m a » 0(1 MeV - 10 GeV) 
/ a « 0(1 - 100 TeV) 


B or F 


No 


d = Q at M* 


Quasi-stable dark matter with: 
m DM « O(10 TeV) 
r DM « O(10 26 sec) 



Table 1: The symmetry structure necessary to realize the scenario presented in this paper. A 
light axion-like state emerges from spontaneous breaking of an R or PQ (Peccei-Quinn) symmetry, 
which then mixes with the Higgs sector of the standard model. Composite states in the strong 
sector are quasi-stable because of a B ( "baryon number" ) or F ( "flavor" ) symmetry. The R or PQ 
symmetry is explicitly broken by supergravity effects or by dimension five operators suppressed 
by Mi ~ O(10 9 - 10 18 GeV). This gives a sufficiently large mass to the axion-like state. Explicit 
breaking of B or F is due to dimension six operators suppressed by M* « O(10 16 - 10 18 GeV), 
leading to dark matter decay through the light states. 

typically involves (long) cascade chains. Second, the mass of dark matter is rather large, 
7Bdm ~ O(10 TeV), since it arises as a quasi-stable state in the supersymmetry breaking 
sector. These lead to a rather broad structure in the electron plus positron spectrum in the 
sub- TeV region after propagation to the earth. Remarkably, we find that such a structure 
beautifully reproduces the spectra recently reported by the FERMI and H.E.S.S. [221 E3] 
experiments. 

Collider signals of light axion-like states: The couplings between the Higgs and axion-like 
states imply that the axion-like states also couple to the standard model gauge fields at 
the loop level. Such couplings may also arise from contributions from the supersymmetry 
breaking sector. This leads to the possibility of producing the light states at the LHC, which 
subsequently decay into standard model fields through mixings with the Higgs fields [24|. 
The final states are most likely leptons, if recent cosmic ray data are explained as described 
above. This may provide a relatively clean signal to discover the light states. 
The couplings of the Higgs to light states also raises the possibility that the Higgs boson 
decays into two axion-like states. If the axion-like state decays mainly into two leptons, then 
this leads to a four lepton final state, whose invariant mass peaks at the Higgs boson mass. 
If the rate is sufficiently large, this leads to a way of seeing the axion-like state (and the 
Higgs boson) at hadron colliders. 

The symmetry structure necessary to produce the signatures described above is summarized 
in Table [TJ While this structure is strongly motivated by the low energy supersymmetry breaking 
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scenario, all that is actually required is some strong dynamics at ps O(10 - 100 TeV) satisfying the 
properties given in the table. Given that QCD already has more or less all the desired ingredients, 
we expect that the required structure may arise naturally in wide classes of strongly interacting 
gauge theories, including non-supersymmetric theories. (For a non-supersymmetric theory, the 
symmetry leading to a light state must be a PQ symmetry. This implies that the theory must 
have two Higgs doublets.) The dynamical scale of O(10 - 100 TeV) then suggests that this 
sector is related to the weak scale through a loop factor. This feature is automatic in low energy 
supersymmetry breaking theories. 

A remarkable thing is that the first signature may have already been seen in the recent cosmic 
ray electron/positron data. The PAMELA experiment found an unexpected rise in the positron 
fraction in the energy range between 10 and 100 GeV, while the FERMI experiment saw an 
excess of the electron plus positron flux over standard diffuse cosmic ray backgrounds in the sub- 
TeV region. These results suggest a new source of primary electrons and positrons with a broad 
spectrum extending up to a few TeV. We find that these features are very well explained by dark 
matter in our framework: a quasi-stable state with mass of O(10 TeV), cascading into leptons 
through light axion-like states with lifetime of O(10 26 sec). We perform a detailed analysis for 
the cosmic ray data and find that a wide range for the mass is allowed for the axion-like state: 
it can take any value between 2m e ~ 1.0 MeV and 2mf, ~ 8.4 GeV except for a small window 
between 2m p ~ 1.9 GeV and 2m T ~ 3.6 GeV. 

These data, therefore, point to the setup of Tabled] as a new paradigm for dark matter, which 
can be beautifully realized in the framework of low energy supersymmetry breaking. Since the 
precise structure of the supersymmetry breaking sector is highly model dependent, one might 
worry that the signatures considered here may depend on many details of the supersymmetry 
breaking sector, which leads to large uncertainties. This is, however, not the case. As emphasized 
above, the existence of the signatures depends only on basic symmetry properties of the super- 
symmetry breaking sector, and their characteristics are determined only by a few parameters such 
as the mass of dark matter and the mass and decay constant of the axion-like state. While we will 
provide an illustrative model as a proof-of-concept, many details of the model are unimportant 
for the signatures. Of course, the flip side of this is that we cannot probe the detailed structure 
of the supersymmetry breaking sector solely by studying these signatures. We may, however, still 
explore some features by carefully studying cosmic ray spectra. 

The signatures described here are complementary to the information we can obtain in other 
methods. In the framework of low energy supersymmetry breaking, the LHC will be able to 
measure some of the superparticle masses. This, however, may not determine, e.g., the scale of 
supersymmetry breaking, since the most general supersymmetry breaking sector provides little 
definite prediction on superparticle masses, as recently elucidated in Ref. [25] • The existence of 
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the very light gravitino can give specific signals, for example, those in Ref. [26j. The signatures 
considered here can add even more handles. In addition to indicating the specific symmetry 
structure of Table [TJ different final states for the axion-like state decay may also be discriminated, 
e.g., by future measurements of the diffuse 7-ray flux at FERMI. These will provide valuable 
information in exploring the structure of the supersymmetry breaking sector. 

The organization of the paper is as follows. In the next section, we describe our supersym- 
metric setup in detail. We explain that quasi-stable states with the desired lifetimes and light 
axion-like states with the desired masses can naturally arise. We also discuss constraints on 
axion-like states, and find that a wide range for the masses and decay constants are experimen- 
tally viable and can lead to leptonic decays. In Section [3j we present an example model that 
illustrates some of these general points. Dark matter is a stable "meson" state in the hidden 
sector that decays into R axions with a lifetime of O(10 26 sec). In Section HI we perform a de- 
tailed analysis of the recent cosmic ray data, and find that the results of PAMELA, FERMI, and 
H.E.S.S. are very well explained. We present a general analysis in the case where dark matter 
decays into e + e~, 7r + 7r~7r°, or r + r~ either directly or through 1-step or 2-step cascades. 

Implications for future diffuse 7-ray measurements are also discussed. In Section [51 we briefly 
discuss collider signatures associated with the existence of light states. Finally, discussion and 
conclusions are given in Section [61 where we mention related alternative scenarios. 

2 Framework 

We consider a supersymmetry breaking sector which consists of fields and interactions character- 
ized by a scale 

A « O(10 - 100 TeV). (3) 

The actual spectrum of this sector could span an order of magnitude or so due to its nontrivial 
structure. The scale of Eq. Q is supposed to arise dynamically through some strong gauge 
interactions in order to explain why the supersymmetry breaking scale, and thus the weak scale, 
is hierarchically smaller than the Planck scale |27j . We assume the sector contains fields charged 
under the standard model gauge group which directly feel supersymmetry breaking. Standard 
model gauge loops then generate gaugino and scalar masses in the super symmetric standard 
model (SSM) sector through gauge mediation [8, 9J. 

The supersymmetry breaking sector may also directly interact with the Higgs fields in the 
superpotential. A possible form for these interactions is 

W = X u H u O u + X d H d O d} (4) 

or 

W = XNH u H d + NO N} (5) 
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Figure 1: A schematic depiction of the setup. 



where H u ^ and N are two Higgs doublet and singlet chiral superfields, respectively, and O u ^,n 
represent operators in the supersymmetry breaking sector. These interactions can generate the \x 
term, and thus lead to realistic electroweak symmetry breaking [10j[TT]. A schematic picture for 
the current setup can be seen in Figure [TJ 

We assume that the entire system of Figure [lj including both the supersymmetry breaking 
and SSM sectors (as well as possible direct Higgs interactions), respects some approximate global 
symmetry under which both the supersymmetry breaking and SSM sector fields are charged. This 
implies that the two sectors must have an interaction that can transmit charges of the global 
symmetry from one sector to the other. We find that an R symmetry is a natural candidate for 
such a symmetry, since interactions generating gaugino masses must always transmit R charges. 
Another simple candidate is a PQ symmetry. In the presence of direct Higgs couplings to the 
supersymmetry breaking sector as in Eqs. ([HE]), charges of the PQ symmetry can be transmitted. 
Note that R and PQ symmetries are global symmetries of the SSM sector only in the limit of 
vanishing /i term and gaugino (and holomorphic supersymmetry breaking) masses. To generate 
the \i term and gaugino masses, therefore, these symmetries must be spontaneously broken by 
the dynamics of the supersymmetry breaking sector. 

The fundamental supersymmetry breaking scale of this class of theories is of order the dy- 
namical scale 

VF^A, (6) 

yielding a light gravitino, m 3 / 2 ~ A 2 /Mpi rs 0(0.1 - 10 eV). This solves all the cosmological 
problems in supersymmetric theories unless the supersymmetry breaking sector introduces its 
own problemsp which we assume not to be the case. The thermal history of the universe is 
normal up to a very high temperature T ^> A, such that the relic abundances of the quasi-stable 
states are determined by a standard thermal freezeout calculation. 

1 Possible problems include the system being trapped in the wrong vacuum or the appearance of stable charged 
or colored particles. 
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2.1 Quasi-stable states 

In addition to the R and PQ symmetries described above, the supersymmetry breaking and SSM 
sectors can have additional independent global symmetries. For example, the SSM sector has 
a baryon number symmetry, at least if R parity is conserved (which we need not assume here). 
Similarly, the supersymmetry breaking sector may possess accidental global symmetries, and since 
this sector is assumed to feel strong gauge interactions, natural possibilities are "baryon number" 
and "flavor" symmetries. If the global symmetry is not spontaneously broken by the dynamics 
of this sector, the lightest state charged under that symmetry is stable. This state can then be 
dark matter, with an abundance determined by its annihilation cross section. 

Since dark matter is the lightest state charged under a global symmetry, it sits at the lowest 
edge of the spectrum of the corresponding charged states. It is therefore natural to expect that 
the dark matter mass is in the range 

m DM ~ O(10 TeV). (7) 

In fact, we will see later that a mass of this size reproduces well the observed electron/positron 
spectrum at PAMELA, FERMI, and H.E.S.S. experiments. The annihilation cross section of dark 
matter is naturally of order 

M ~ ^4-- (8) 

where k represents typical couplings between states in the strong sector, and we have assumed 
two-body annihilation. For n pa 3, this gives the cross section needed to reproduce the observed 
dark matter abundance, Odm — 0.2. Such a value for k is quite natural for couplings between 
hadronic states in a strongly interacting sector. 

The global symmetry ensuring the stability of dark matter is expected to be an accidental 
symmetry at low energies. This implies that the symmetry is not respected by physics at some 
higher energy M*, such as the unification or Planck scale, so that the effective theory at the 
scale A contains higher dimension operators suppressed by powers of M* that do not respect the 
global symmetry. This situation is precisely analogous to baryon number in the standard model 
embedded in grand unified theories (with M* identified with the unification scale), or strangeness 
in QCD (with M* identified with the weak scale). Dark matter can then decay with cosmologically 
observable timescales, depending on the dimension of the leading symmetry-violating operators. 

Consistent with gauge coupling unification, we take M* to be around the unification or Planck 
scale, M* pa O(10 16 - 10 18 GeV). Then, to have observable signatures in cosmic rays, the dimen- 
sion of relevant symmetry violating operators must be six, as illustrated in Ref. [20] in the case 
where m^M ~ O(TeV) and M* pa O(10 16 GeV). With dimension six decay operators, we find 

o M t n m25 ( M * \V 10TeV \ 5 

tdm « 8vr^ ~ 2 x 10 25 sec — — — , 9 

m^ M VIO 1 ' GeV/ V m DM / 
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for two-body decays, and 

., Ml 97 ( M* \ 4 /10TeV\ 5 , , 

r DM re 128tt 3 ^- ~ 3 x 10 27 sec * — ) , 10 

m^ M VIO 1 ' GeV/ V m DM / 

for three-body decays. The required lifetime to fit the PAMELA, FERMI, and H.E.S.S. data 
through dark matter decay is of order 10 26 sec, which is consistent with the value of M* taken 
here. 



2.2 Decay of quasi-stable states 

There are several possible ways for the decay of dark matter to be caused by dimension six 
operators. One simply needs small breaking of the global symmetries protecting dark matter, and 
generically decays will proceed via some kinematically-allowed but symmetry-violating channel. 
In particular, dark matter can decay into light axion-like states and/or the gravitino. As long 
as the dimension six operators are suppressed by the unification or Planck scale, they can even 
explicitly violate the R or PQ symmetry that gives rise to the axion-like state, since such explicit 
breaking gives very small masses compared to the contributions considered in Section 12.31 

Suppose that interactions in the supersymmetry breaking sector are asymptotically free, so 
that the dimensions of the operators are determined by the canonical dimensions of the elementary 
fields. Suppose also, for illustrative purposes, that the supersymmetry breaking sector contains a 
gauge group SU(N C ) (N c > 3) with N f flavor of "quark" fields Q' + Qi (i, % = 1, • ■ ■ , N f ), where Q i 
and Qi are chiral superfields in the fundamental and anti-fundamental representations of SU(N C ), 
respectively. We also assume that the quark fields have generic masses of order A (or somewhat 
smaller) due to some (i?-violating) dynamics in the supersymmetry breaking sector. 

The SU (N c ) gauge group is supposed to confine at a scale ~ A. This then leads to composite 
"meson" fields M| ~ Q l Qj, and also "baryon" and "antibaryon" fields B ~ Q Nc and B ~ Q Nc for 
Nf > N c . Because of the nonzero quark masses, these fields have masses of order A (or somewhat 
smaller). Moreover, in the limit that nonrenormalizable operators vanish, the states Mj (i ^ j), 
B and B can easily be stable (except for B and B for N c = 3). These states are therefore good 
candidates for dark matter. 

The lifetimes of the quasi-stable states described above are controlled by the form of higher 
dimension operators. Let us begin with the case where dark matter is identified with baryons 
(and antibaryons) . The decay of dark matter can then occur through baryon number violating 
operators in the superpotential. These operators are dimension six if N c = 5% 

W ~ -^QQQQQ + -LqQQQQ, (11) 

2 In the following equations, any one of the operators is sufficient to cause dark matter decay. In particular, 
physics at M* need not respect a Q <-> Q symmetry. 
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in which case the lifetime of dark matter can be of order 10 26 sec as needed to reproduce the 
electron/positron data§| 

An alternative possibility is that dark matter decays through holomorphic terms in the Kahler 
potential K, which can also be viewed as superpotential terms suppressed by an extra power of 
Mpi after performing a Kahler transformation. The required lifetime is obtained for iV c = 3: 

K~±QQQ + ± r QQQ + h. c . =► W ~-fj—QQQ + -f—QQQ, (12) 

where we have used the fact that the superpotential has a constant term of order A 2 Mpi to cancel 
the cosmological constant, and assumed that the renormalizable superpotential term W ~ QQQ 
is absent, perhaps because of an R symmetry. Note that in both of the cases above, the number 
of "colors" N c needs to be chosen appropriately to have dimension six dark matter decay. 

We now consider the case where the meson states M| (i ^ j) are dark matter. These states 
exist even for Nf < N c , and the longevity of their lifetime could be ensured by a vector-like 
U(l) Nf symmetry that may exist in the Lagrangian at the renormalizable level. The decay of 
these states can be caused by dimension six operators in the Kahler potential: 

K ~ ^Q\Q j QlQ l + j^Q r %Q fk Qi+ j^Q\Q 3 Q fk Qi. (13) 

Possible lower dimension operators in the superpotential W ~ (l/M^)Q l QjQ k Qi can easily be 
absent, for example by imposing an R symmetry. An attractive feature of this possibility is that 
the lifetime does not depend on the number of colors N c , and that the number of flavors Nf need 
not be equal to or larger than N c . This setup, therefore, can naturally be accommodated in a 
wide variety of gauge theories. In Section [3] we present an explicit model realizing this possibility, 
where the quasi-stable mesons decay via Eq. (fl3"j) into R axions. 

Although we have not considered them in the discussion above, in general the supersym- 
metry breaking sector also contains states not charged under SU(N C ). Dark matter decay 
may then occur through these states. For example, mesons may decay through operators of 
the form K ~ QlQi & $ / M% + Q^Qj&Q/M^, where $ represents a generic field in the su- 
persymmetry breaking sector. The existence of singlet fields can also change the requirement 
on N c for the baryon dark matter case. As long as the lowest symmetry breaking opera- 
tors are dimension six, however, the existence of these operators do not affect the basic argu- 



3 A long-lived hidden sector baryon with the SU (5) gauge group was considered in Ref. [28] in the context that 
the quarks Q l are also charged under the standard model gauge group. This model, however, does not preserve 
the success of perturbative gauge coupling unification, since the extra matter content charged under the standard 
model gauge group is too large, see e.g. [25]. Here we consider that the quarks Q l are not charged under the 
standard model gauge group. 



10 



mentjj Below, we assume for simplicity that operators containing only quarks dominate the 
decay. Whether this is the case or not is determined by physics at the scale M*. This also 
implies that direct interactions between the supersymmetry breaking and SSM sectors, such as 
K ~ QIQ3$I su $ S sm/M2 + Q^Qj$I su $ S sm/M%, are relatively suppressed. Here, $ S sm represents 
SSM fields. We will comment on the case where these interactions are relevant in Section [6j 

2.3 Light axion-like states 

Spontaneous breaking of an approximate R or PQ symmetry in the supersymmetry breaking 
sector leads to a pseudo Nambu-Goldstone boson. The phenomenology associated with this 
particle is determined largely by its mass m a and decay constant f a . Suppose that relevant 
interactions in the supersymmetry breaking sector obey naive dimensional analysis for a generic 
strongly coupled theory [30J. In this case, we find that the decay constant is naturally somewhat 
(« 47r) smaller than the dynamical scale A « O(10 - 100 TeV): 

/ « 0(1-10 TeV). (14) 

According to naive dimensional analysis, the generic size for expectation values of the lowest 
and highest components of a chiral superfield $ are given by ($) ~ A/Att and (F$) ~ A 2 /47r, 
respectively, and a generic coupling constant has the size k « Atc. (Here we have ignored a 
possible N c or Nf factor associated with the multiplicity of fields, but this does not affect the 
basic discussion.) This implies that a generic supersymmetric mass and supersymmetry breaking 
mass-squared splitting in the supersymmetry breaking sector are of order M mess « « A and 
F mess ~ K (F<!>) ~ A 2 , respectively, so that the gaugino and scalar masses generated in the SSM 
sector are of order (g 2 /167r 2 )F mess /M mess « (g 2 /167r 2 )A, which is consistent with Eq. (j3J). On the 
other hand, the decay constant of an axion-like state scales as f a ~ ($) ~ A/47T, giving Eq. (Tl4|) . 
Note that the suppression of f a over A here is precisely analogous to the fact that in QCD the 
pion decay constant, f w ~ 0(100 MeV), is an order of magnitude smaller than the characteristic 
QCD scale, i.e. the rho meson mass m p « 0(1 GeV)§| 

The mass of a light axion-like state is determined by the size of explicit symmetry breaking. 
Let us first consider the state associated with an R symmetry — an R axion. An interesting 
feature of an R axion is that it has an irreducible contribution to its mass from supergravity, 

4 With singlets, some of the accidental symmetries are not as "automatic" as the case without singlets. For 
example, the low energy U(l) Nf flavor symmetry does not exist if the quarks couple to two or more singlets with 
arbitrary Yukawa couplings. 

5 In general, it may not be true that all the couplings in the supersymmetry breaking sector are strong. In this 
case, the decay constant need not obey Eq. ([M]) : for example, it can be easily of O(100 TeV). In particular, if the 
messenger fields (fields in the supersymmetry breaking sector charged under the standard model gauge group) feel 
i?-breaking effects though perturbative interactions, then f a for an R axion needs to be of O(10 TeV) or larger to 
generate sufficiently large gaugino masses. 
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whose size can be determined by the strength of fundamental supersymmetry breaking F. This 
is because to cancel the cosmological constant, the superpotential must have a constant piece 
(W) Rj FMpi ps 0(A 2 M P1 ), which is necessarily R- violating. Since (W) ^> A 3 , the constant piece 
must come from a sector other than the supersymmetry breaking sector, implying that it appears 
as explicit breaking from the perspective of the supersymmetry breaking sector. This provides 
the following contribution to the R axion mass [T5] 



where (<&) is a generic vacuum expectation value in the supersymmetry breaking sector, and 
we have used (<&) ~ A/4tt, f a ps A/Att, and F ~ A 2 /Att. The uncertainty of the estimate, 
however, is very large, so that we can easily imagine m a in the range 0(1 - 100 MeV). Unless 
the supersymmetry breaking sector contains another explicit breaking of the R symmetry, the 
R axion mass is given by Eq. (|T5|) . 

We now consider the case of a PQ symmetry, or of additional explicit breaking of an R 
symmetry. If the symmetry is violated by dimension five operators in the supersymmetry breaking 
sector, then we expect 




where we have introduced the scale Mj for physics causing the explicit breaking. We imagine Mj 
to take a value between an intermediate scale and the gravitational scale, Mj ps O(10 9 - 10 18 GeV), 
giving m a ps 0(1 MeV - 10 GeV). While smaller masses are possible, they are in conflict with 
astrophysical measurements as we will see later. The origin of the intermediate scale might 
be associated, for example, with the Peccei-Quinn scale for the QCD axion, the scale of the 
constant term in the superpotential, or the B — L breaking scale generating right-handed neutrino 
masses^ Alternatively, the explicit breaking may be due to a tiny dimensionless coupling in 
the supersymmetry breaking sector. Note that a possible contribution from the QCD anomaly, 
m a ^ m n (f n /f a ) ps 0(1 keV)(100 TeV/A), is small. 

2.4 Couplings of axion-like states 

The couplings of the axion-like states to the SSM sector are completely determined by the sym- 
metry structure. We are considering the case in which the Higgs bilinear h u hd is charged under 
the symmetry that leads to the light state, where h U) d are the lowest components of H u ^- This is 

6 If we use f a k, A instead of A/4tt, we obtain m 2 a w O(100 MeV) 2 (A/100 TeV) 3 (10 17 GeV/Mj). This gives 
m a w 0(1 MeV - 1 GeV) for Mj w M*. 




Mpi/2 4ttM p1 



F($) 3 A 3 



O(10 MeV) 2 ( , 

v ; VlOOTeV/ 



(15) 
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almost always true for an R symmetry (unless direct couplings of the Higgs to the supersymmetry 
breaking sector force vanishing charges for H Ujd ), and is by definition true for a PQ symmetry. In 
both cases, we can take a field basis in which the axion-like state a is mixed into the h„ d fields: 



u,d 

„■ COB 2 /3 „ 

K = v u e^u a , (17) 
h d = v d e^ a , (18) 

where v Uyd = (h u> d) and tan/3 = (h u )/(hd). The distribution of a inside h Uyd is determined by 
the condition that a is orthogonal to the mode absorbed by the Z boson, and we have arbitrary 
chosen the 0(1) normalization for f a . The expressions of Eqs. ( !T7l [TBI) completely determine 
the leading-order couplings of a to the standard model quarks and leptons. For example, the 
couplings to the up quark, down quark, and electron relevant for a decay are given by 

.m u cos 2 8 - ,m d sm 2 3 - m e sm 2 3 - 

C=-t u / -^ a^ ul5 ^ u -t / am dlb m d -i M afr e75 fr e , (19) 

V 2/ a V Ija V ^ Ja 

where m u ^ e are the up quark, down quark, and electron masses. The couplings to heavier 
generation fermions are similar. 

The couplings of a to the Higgs boson h can be obtained by replacing v Uyd as 

cos a sin a 

v u -> v u + —/^h, v d ->v d - —f^h, (20) 

in Eqs. < \I7\ [T8|) and plugging the resulting expressions into the Higgs kinetic terms C = \d^h u \ 2 + 
\d^h d \ 2 . Here, a is the Higgs mixing angle. This leads to 

C = ^- 2 h(d fl a) 2 + ^- 2 h 2 (d,a) 2 , (21) 



where v — Jv 2 + v d , c\ — sin 3 cos /5(cos 3 /3 cos a — sin 3 3 sin a) , and c 2 = cos 4 /5 cos 2 a + sin 4 /3 sin 2 a;. 
In the decoupling limit, a ~ 3 — ir/2, this equation reduces to 

fsin 2 2/?, / „ , 9 sin 2 2/? 9/ „ , 9 , , 

£ = W + -y^r * (V) ■ ( 22 ) 

The first term is responsible for Higgs decay into two axion-like states. 

There are also couplings between a and the standard model gauge bosons. Their precise values 
depend on the matter content of the entire theory. At a given energy scale E, effective direct 
interactions between a and the gauge bosons are given by 

C = Y ^ aF^ F AfJjU , (23) 
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where A runs over the gauge groups accessible at the scale E, and g& and ca are the gauge coupling 
and a coefficient of order unity. The coefficient ca encodes the contributions from physics above E, 
and is determined by the symmetry properties of the fields integrated out to obtain the effective 
theory At the electroweak scale or below, we generically expect ca ^ for color, A = SU(3)c, 
and electromagnetism, A = Z7(1)em- The interaction with gluons (A = SU(3)c) is responsible 
for direct production of a at hadron colliders. 

2.5 Constraints on axion-like states 

The axion-like states for the R and PQ symmetries considered here couple to the standard model 
fields as the DFSZ axion [31j . except for possible differences in the numerical coefficients ca- 
Their masses, however, are much heavier because of explicit symmetry breaking, so that the 
experimental constraints on them are quite different from those on the DFSZ axion. Here we 
summarize the constraints on m a and f a for the region relevant to our discussions. For previous 
related analyses, see Refs. [TSl [32] • 

The constraints on an axion-like state a depend strongly on its decay mode. The decay width 
of a into two fermions is given by 

Fa-// ) = J-L-!-l ll-—f) > 24 
16tt / a 2 \ mi J 

where n/ = 1 and 3 for leptons and quarks, respectively, and c/ = sin 2 /? for / = e, fi,r,d, s,b 
while Cf = cos 2 /3 for / = u,c, t. The decay width into two photons is 

2 4 ^ 
C 6 777 

r < a ^-3S(kv7f' (25) 

where c 7 represents the ca coefficient for £7(1) EM at energies below m e , and e is the electromagnetic 
gauge coupling. 

For m a < 2m e , a decays mainly into two photons with the decay width of Eq. (|25|) . giving 

A possible decay of an R axion into two gravitinos has the width T(a — > GO) w (l/167r)m^ 2 m a / f 2 , 
and is thus negligible unless m a ^ (47r/a)m 3 /2- For 2m e ^ m a ^ 2m M , a decays dominantly into 
e + e~ with the decay width of Eq. (|24"I) . giving 



sm p V m, 



-1/2 
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For 2m^ < a ^ 800 MeV, the a — > mode dominates with 

4 1 /300 MeV\ (fa Y(, ^IY 1/2 , oa . 

In this region, the a — > 7r7r mode is suppressed by CP invariance, and a — > 7T7T7r has the width 
of order (l/128ir 3 )(m^m a / f^f^ ), which is much smaller than T(a — > fi + fi~). For m a <: 800 MeV, 
the a — > p*7r — > 7T7T7r and a — > r]7T7r modes become important, but the final states still contain a 
significant fraction of leptons from charged pion decay, unless m a ^ 2 GeV where nucleon modes 
start dominating. For 2m T < m a < 2m;,, a will decay dominantly into taus with 

cw^ = 6.3 x 10- m > (-^V (l - ^V" (29) 

sin p V m a / \10 TeV J \ mi J 

The branching ratio into cc is ~ ?>m\jm 2 T tan 4 /3, which is highly suppressed for tan/5 ^ 2. For 
m a > 2m the a — > bb mode dominates. 

Rare decays of mesons provide strong constraints on axion-like states. In particular, the 
K + — > 7r + a process gives significant constraints on the region of interest. The theoretical estimate 
for the branching ratio is Br(K + -> vr+a) Z 1.1 x 1(T 8 (10 TeV// a ) 2 [33J. For m a < 2m M , the a 
decay length is large enough that a appears as an invisible particle in K decay experiments. The 
experimental bound on the branching ratio is then ~Bt(K + — ► n + a) ^ 7.3 x 10~ n [31]. While the 
theoretical estimate has large uncertainties, this gives the rough bound 

f a > 100 TeV for m a < 2m M . (30) 

For 2m ^ < m a < tuk — m n , a decays quickly into so that the relevant experimental data 

is Br(K + — > 7r + / u + /i _ ) ~ 1 x 10~ 7 [35], which is consistent with standard model expectations. 
Considering that the dimuon invariant mass would be peaked at m a for a decay, Br(K + — ► n + a) 
should be somewhat smaller than this number, giving a conservative bound of 

f a <; a few TeV for 2m^ < m a < m K — m n . (31) 

Radiative decay of T also provides constraints, but the bounds are typically f a ^ O(TeV) (for 
a — ► [36]) or weaker (for a — ► t + t~ [37]). and do not significantly constrain the parameter 

region considered here. 

There are constraints from beam-dump experiments. For example, the experiment of Ref. [38] 
excludes f a ^ (10 - 100) TeV for 2m e < m a < 2m, M . None of these experiments, however, gives 
as strong bounds as the ones from kaon decay given above, except for possible small islands in 
parameter space. Constraints from reactor experiments are also similar. They are not as strong as 
those from kaon decay, except that the experiment of Ref. [39] excludes a region of f a somewhat 
above 100 TeV for 2m e < m a < 10 MeV. 
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2m e rriK — 'm-K ~ 800MeV 2m p 2m T 2m\, 



a ~ *■ 77 e + e 7r + 7r 7r° hadronic r + r 66 

Figure 2: A schematic picture for the constraints on the m a -f a plane, with the shaded region 
corresponding to the excluded region. Note that the actual limits on /„ have 0(1) uncertainties, 
as explained in the text. The dominant a decay mode for a given value of m a is also depicted. 



Astrophysics provides strong constraints on very light axion-like states [40J. A combination 
of the bounds from the dynamics of the sun, white dwarfs, and horizontal branch stars excludes 
m a ^ 300 keV for the relevant region of f a ~ 0(1 - 100 TeV). Supernova 1987A could also 
provide potentially strong constraints. If an a produced in a supernovae were to freely escape, 
then its mass would be excluded up to ss 0(1 GeV) with the a production rates corresponding 
to the f a values considered here/| This is, however, not the case. For / w 0(1 - 100 TeV), the 
produced a is either trapped inside the supernovae or decays quickly, so that it does not carry 
away significant energy. In particular, for m a > 2m M , a immediately decays into muons, which are 
then thermalized quickly. Supernova 1987A, therefore, does not strongly constrain the parameter 
space considered here, except that the region m a < 2m e is excluded by nonobservation of 7 rays 
from a decays [4"I] . 

A schematic diagram summarizing the bounds on m a and /„ is depicted in Figure [2j together 
with the dominant a decay modes. For f a w 0(1 - 100 TeV), the most natural region avoiding 
all the constraints are m a > 277^, although m a < 2m e is still allowed for f a w O(100 TeV). For 
the range m a rs 0(1 MeV - 10 GeV) considered in Eqs. (fT5j [T6|) . a natural a decay mode will 
be either r + r~, 7r + 7r~7r°, or hadronic (or e + e~ for f a m O(100 TeV)). As we will see in 

Section HI a decay into r + r _ , or 7r + 7r~7r° (or e + e~) can provide a spectacular explanation 

for the recently observed electron/positron signals in cosmic rays. 



7 We thank Savas Dimopoulos for discussion on this point. 
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3 Illustration 



In this section, we present an example model which illustrates some of the general points discussed 
above. The model has quasi-stable meson fields carrying nontrivial flavor, which we identify 
with dark matter. We find that this dark matter, in fact, decays dominantly into R axions (an 
R axion and a gravitino if the lightest meson field is a fermion) through dimension six operators. 
The lifetime of dark matter can therefore be naturally of O(10 26 sec). The resulting R axion 
decays into standard model particles as shown in Section I2.4[ leading to astrophysical signatures 
discussed below in Section HJ We also show that under reasonable assumptions, the correct 
thermal abundance for dark matter can be obtained through a standard freezeout calculation. 

3.1 Setup in supersymmetric QCD 

Ideally, one would want to construct a complete model where supersymmetry (and R symmetry) 
is dynamically broken and successfully mediated to the SSM sector. It is, however, notoriously 
difficult to realize such explicit constructions |42j. Here, we will use the example of supersym- 
metric QCD to show that quasi-stable dark matter candidates can arise from strong dynamics 
with spontaneous i?-symmetry breaking. In Section 13.44 we will show that this strong sector 
can be consistently coupled to supersymmetry breaking dynamics and to messenger fields, and 
we will comment on the parametric difference between the dark matter mass and the scale of 
supersymmetry breaking. 

Our starting point here is supersymmetric QCD where the number of flavors Nf is less than 
the number of colors N c . In the ultraviolet, the matter content consists of quark fields Q l and Qi 
(i,T — 1, • • • , Nf). We couple these fields to a singlet field S in the superpotential: 

W = XlSQ'Qj. (32) 

With only one singlet, these couplings can always be diagonalized by rotations of the quark fields, 
\\ = Xi5j, and we can then easily see that this theory has an accidental U(l) Nf flavor symmetry, 
where each U(l) factor corresponds to a vector-like rotation of the quark fields of a given flavor. 
The theory also possesses a non-anomalous U(1)r symmetry with 

N 2N 

R(Q) = R(Q) = i - jf , R(S) = (33) 

To realize our setup, the R symmetry must be spontaneously broken. Here we treat its effect 
through the background expectation value of the S field, which is sufficient to understand the 
properties of dark matter. We will eventually allow S to be a dynamical field in Section 13.41 
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Below the dynamical scale of QCD, which we take to be close to the scale charactering the 
entire supersymmetry breaking sector, § the appropriate degrees of freedom to describe the dynam- 
ics are composite meson fields M| ~ Q l Qj interacting via a non-perturbative superpotential [13]. 
Together with the couplings of Eq. fl32l) . this leads to the effective superpotential 

/ \Nf \ N c -N f 

W = X t ASM t + X M (N c -N f )A 3 {^-^\ , (34) 

where A is the dynamical scale, and Mj = M\ (i = i) are the diagonal mesons. The coefficient 
Am is an unknown factor coming from canonically normalizing the meson fields, and a similar 
factor for the first term is absorbed into the definition of Aj. Using naive dimensional analysis, 
the size of these coefficients are X M « O^Tr)^ 27 ^-^/^"^) and Aj/Ai « 0(1/4tt). Note that 
there are no low energy baryon fields when Nf < N c . 

Setting (S) ^ but (Fs) = 0, the superpotential of Eq. ( 1341) has a stable supersymmetry- 
preserving, R- violating minimum. Defining = Xi(S) and detm = fli^j; the minimum of the 
potential is at 

(M 4 ) = ^, (M;) = (i^j), (35) 

TTli 

where for convenience we have defined olm = Aj^(det m/X^/A N f ) 1//Ar<: . According to naive dimen- 
sional analysis, mj ps 0(X{S) / An) and ~ 0((X(S) /A) N f/ Nc /I6ir 2 ), where A represents the size 
of the original couplings in Eq. fl32l) . Since (M|) = for i ^ J, the U(l) Nf flavor symmetry is 
unbroken, making the Mj mesons with i ^ J (for now, absolutely) stable. Therefore, the lightest 
components of these fields, either scalars or fermions, are dark matter candidates. 

The diagonal mesons Mj are unstable. Since the R symmetry is assumed to be spontaneously 
broken, the theory has a light R axion with the decay constant f a of order A. The couplings of the 
mesons to the R axion are determined by symmetries, and by doing appropriate field redefinitions, 
these couplings can be read off from the meson kinetic termsj^l For the scalar components, the 
kinetic term is C = \d^Ma + £cAf((Mj) + Mjj)(<9 M a)// a | 2 , where Ma = Mj — (Mj) and cm is an 
0(1) coefficient that depends on R charges of the fields. Assuming real (Mj), this gives a coupling 
of Re(Mjj) to two R axions, and coupling between Re(Mjj), Im(Mjj), and an R axion. Therefore, 
Re(Mjj) decays promptly into two R axions, while Im(Mjj) decays (somewhat less) promptly into 
three R axions via an off-shell Re(Mjj). After supersymmetry breaking, the fermionic component 
of Mj-j has a coupling to the corresponding scalar component and a gravitino, suppressed by 
powers of A. This allows the Ma fermion to decay promptly into a gravitino and an R axion. 



8 We discuss a potential difference between the dynamical scale A appearing here and the scale that breaks 
supersymmetry in Section 13.41 

9 Once the S field becomes dynamical, S will mix with the meson fields. This can introduce extra contributions 
to the couplings, especially if there is no hierarchy between the scales involved. It, however, does not affect the 
basic conclusion here. 
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3.2 Thermal relic abundance 



Expanding around the minimum of the potential from Eq. f[35|) . the mass terms for the mesons 
are 

W = 2^ £ (^M,,^ + __M,,M^ j , (36) 

where we have used the notation My = (Mj — (M|))|j=y. The stable mesons (My with i 7^ j) 
do not have mixing terms because of the unbroken U(l) N f symmetry, while the unstable mesons 
(Mjj) have a mixing term that is suppressed by l/(N c — Nf). For simplicity of discussion, we will 
focus on the large N c limit where there is a simple relationship between the stable and unstable 
meson masses: 



mass(Mjj) ~ y mass(Mjj) mass(M_y). (37) 
Up to l/(iV c — Nf) corrections, the leading interaction term for the mesons is 

W = =- Y MuMikMki \i fc + — =- Y MijMj k M kl Mu \ . 38 

3a M~t A 3 4a^^J J J A 5 v ; 

The superpotential coupling in Eq. ( l38l) allows for various annihilation diagrams for stable 
My mesons. By the mass relation in Eq. (|37|) . either the ii-type or jj -type mesons will be lighter 
than the zj-type mesons. Assuming rrii < rnj, the annihilation cross section for MyM^ — > MaMa 
scales like 

1/t 4 1 rriimj 1 mfrrij 

(av) ~ — , m DM t— , k ~ — t~o • (39) 

87rm£, M «m A a 2 M A 6 

Using naive dimensional analysis, these quantities can be estimated as 



m D M ~ A 



Nf IN c 

(a) • K ~ 4 n^r) • (40) 



so that the required annihilation rate can naturally be obtained. For example, it is easy to obtain 
wdm ~ O(10 TeV) and k a factor of a few, which leads to (av) ~ (l/87r)(l/TeV 2 ), and thus 
flj)M — 0.2. Note that here we have not included possible multiplicity factors in naive dimensional 
analysis, which will somewhat decrease the value of k. We therefore do not prefer very large N c in 
practice. Note also that, strictly speaking, our analysis is not theoretically very well under control 
for \(S) A, although we expect that the basic dynamics is still as presented for X(S) ~ A. 

There is subtlety when annihilation occurs into states of comparable mass. If ~ rrij, then 
the ij-, ii-, jj -type are nearly degenerate for relatively large N c . In this case, one must use the 
methods of Ref. [2] to properly calculate the thermally averaged annihilation rate. One can, 
however, still obtain the right thermal relic abundance with somewhat stronger couplings. If 
Nf is large, then there is a potential for multiple dark matter components having comparable 
abundances. The freezeout calculation can then be affected by various co-annihilation channels. 
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3.3 Dark matter decays 

The U(l) Nf flavor symmetry which ensures the stability of dark matter arises as an accidental 
symmetry of the renormalizable interactions of Eq. (|32|) . As such, it is plausible that this symme- 
try is not respected by physics at the unification or gravitational scale. Suppose that the leading 
operators encoding the high energy physics are the Kahler potential terms 

K = TpVilQtQ'QW, (41) 

with arbitrary flavor structures for rf^ in the basis where the interactions of Eq. (1321) are diagonal: 
Af = Xi5j. Below the scale A, these operators can then be matched intc0 

K = ^4 M M, (42) 

where cij k are coefficients. Using naive dimensional analysis, we find Cj k ~ 0(l/167r 2 ) for rf-l ~ 
0(1) and \i{S) ~ 0(A). In general, one expects Cj k to have 0(1) OP-violating phases. 

The existence of the terms in Eq. (j42p induces mixings between the diagonal and off-diagonal 
meson states, leading to the decay of dark matter. The decay width has a large suppression from 
the A 2 j Ml factor in Eq. (I42p . so that the lifetime of dark matter is very long. One possible decay 
chain arises through M u —>■ aa, where a is the R axion. Through the mixing, this leads to 

Mij -> aa, (43) 

where we have assumed that dark matter is a scalar^ For fermionic dark matter, one of the 
a's must be replaced by a gravitino G. Another possibility is that dark matter decays through 
the meson self interactions of Eq. (|38[) . Again, through the mixing, this can lead to 

M U M U , (44) 

if it is kinematically allowed, mass(Mjj) > 2mass(Mjj). Here, Ma in the final state can be either 
a scalar or fermion, which subsequently decays into aa or aG, respectively. 

For the scalar dark matter decay in Eq. (|43l) . the lifetime of dark matter can be estimated as 

ft /M 2 \ 2 97 /1/16tt 2 \ 2 / Af, V 4 /10TeV\ 5 , , 

TDM " 8 *^ZW (sp) = 3 x 10 sec ( c ) (lira) (— ) • (45) 

where (M) is a typical meson expectation value, c represents a generic size for the coefficients 
c l L in Eq. (|42|) . and we have used m^M ~ 47r/ a pa A in the last equation. For M* of order the 



10 The same argument as below applies to the operator K = (l/M^)Q^QjQ^ k Qi or (1 / M%)Q\Q j Q* k Q 

Cj l k , even Im(-My 



11 Because of the CP- violating phases in c\ l k , even Im(My ) dominantly decays into two R axions, not three. 
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unification scale, the lifetime is in the range required to produce observable cosmic ray signatures. 
While there are different parametric dependences for tdm if one uses the decay mode in Eq. (jHJ) 
or if one considers fermionic dark matter decay, they all give the same estimate as Eq. (1451) after 
using naive dimensional analysis. 

3.4 Towards a more complete theory 

In the above discussion, we assumed that (S) could be treated as a spurion for the spontaneous 
breaking of the R symmetry. In a complete theory, S would be a propagating degree of freedom 
that is part of a larger supersymmetry breaking sector. Naively coupling 5* to the mesons as in 
Eq. (1341) . the i^-term potential would change the meson vacuum structure. In particular, even 
if S could be stabilized through terms in the Kahler potential, there is generically a mesonic 
runaway direction [45] . 

Therefore, it is important to have a proof-of-concept that a propagating 5* field can not only 
obtain an R- violating vacuum expectation value but also couple to the mesons without introducing 
runaway behavior. In addition, one would like to see that appropriate messenger fields can be 
added to the theory to communicate supersymmetry breaking to the SSM. While a complete 
theory would incorporate these two effects in a completely dynamical setting with all scales set 
by dimensional transmutation, here we consider O'Raifeartaigh-type models in order to treat 
these effects modularly (with a dynamical model mentioned only briefly towards the end). We 
leave further model building to future work, although we emphasize that the precise details of 
supersymmetry breaking are largely irrelevant for the dark matter discussion. 

We will also not be concerned by the specific mass hierarchies needed to realize a realistic 
theory. Typically, the scale of supersymmetry breaking should be closer to 0(100 TeV) to obtain 
a realistic superparticle spectrum, while the mass scale for dark matter suggested by the cosmic 
ray data is closer to O(10 TeV). In a single scale theory with naive dimensional analysis, both 
mass scales are expected to coincide. However, in a complete dynamical theory of supersymme- 
try breaking, there may be additional structures that generate such a "little" hierarchy] 12 ! For 
example, this can be realized if the SU(N C ) theory in Section 13.11 is in a strongly interacting 
conformal window at high energies. If nontrivial dynamics kicks in at the scale O(100 TeV) to 
break supersymmetry, the same dynamics can make the SU(N C ) theory deviate from the fixed 
point, presumably due to decoupling of some degrees of freedom, leading to the dark matter setup 
considered here. The physics of R breaking will be associated with the higher scale dynamics. 

To see how S can be made dynamical, consider adding the following interactions to Eq. f)34p : 



12 This little hierarchy may not be necessary in general if the axion-like state decays into ir + Tr 7r° or t + t (see 
Figured]) and the number of messenger fields is relatively large. 




(46) 
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These interactions could be the low energy description of a more complete supersymmetry break- 
ing sector, and so we consider that [is, m Si and rhs are roughly of order A. Regardless of the 
R charge of the S field (i.e. regardless of Nf and N c ), there is a unique consistent R charge 
assignment for the new fields, with 

2N 2N - 2N 

R(X) = 2, R(S) = -^ R(Y) = 2-^, R(X) = 2 + ^. (47) 

7V/ 7V/ JSf 

With the inclusion of the ms term, the Fg equation of motion simply sets the value of Y and does 
not change the potential for the mesons, so that the global minimum of the meson potential is still 
given by Eq. fl35|) . The Ax coupling forces S and S to obtain vacuum expectation values, and the 
ms and % terms constrain 5* and 5* from running away to infinity. As long as A x /x| > msfhs, 
this theory has a supersymmetry- and i?-breaking global minimum away from the origin, with S 
vacuum expectation value given by 



(S) 



\ 




(48) 



At tree level, the only fields that obtain F-term vacuum expectation values are X, Y, and Y, so 
the off-diagonal mesons only feel supersymmetry breaking at loop level. The vacuum expectation 
value of X is unspecified at tree level, but one expects loop effects to stabilize X, though not 
necessarily at (X) = 0. Since supersymmetry and R are both broken regardless of the value of 
(X), this is an example of a model where supersymmetry and R are spontaneously broken at tree 
levelH 

In principle, one could calculate the full spectrum of the theory at one loop, allowing one to test 
whether the stable dark matter is the fermionic or bosonic component of the off-diagonal mesons. 
In practice, this is a nontrivial exercise given the large amount of mixing between the various 
new degrees of freedom. If the dynamical sector is truly strongly coupled, then loop counting 
would not apply, and the mesons could feel 0(1) supersymmetry breaking effects depending on 
the details of the Kahler potential. Therefore, one would not expect a one-loop calculation to get 
the correct sign for the fermion/boson splitting. 

We have seen that Eq. (132|) can be consistently coupled to the supersymmetry- and /^-breaking 
sector in Eq. (Hoi) . It is then straightforward to couple messenger fields to the supersymmetry 
breaking, using an analogous structure to Eq. (T46I) . Regardless of how X is stabilized, the Y 
field obtains vacuum expectation values in both the lowest and highest components, so Y can 
communicate supersymmetry breaking to the SSM sector through complete SU(5)sm messenger 



13 In the language of Ref. [46], Eq. (j46|) is the "<? = 0" superpotential which preserves a U(1)r symmetry and 
spontaneously breaks a global U(1)a symmetry. The inclusion of Eq. ([32| identifies a linear combination of U (1)r 
and U(1)a as the true U(1)' R symmetry, such that the R symmetry is broken. 
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multiplets. Consider adding the superpotential 

W = X Y YFF + m F FF' + rh F F'F, (49) 

where F and F' are 5s of SU(5i)suj an d F an d F' are 5s, and there is a consistent R charge 
assignment for the messengers. Taking mpmp > Xyms(S), we can ensure that the messengers do 
not develop vacuum expectation values. With this kind of messenger sector, the gaugino masses 
are proportional to F Y instead of Fy [47J. They are, therefore, suppressed compared with the 
scalar masses for XyFy <C mpmp, which may not be fully desired. Moreover, one also needs to 
take rriF — trip to avoid an unwanted large Fayet-Iliopoulos term for U(l) hypercharge. As a 
proof-of-concept, however, we see that it is possible to couple the dark matter, supersymmetry 
breaking, and messenger sectors together in a consistent way. 

If one desires, the above model can be extended in such a way that all the dimensionful scales 
Us, m s , rh s , m F , and rh F are generated dynamically. These scales can be replaced by a single 
chiral superfield T, which obtains a vacuum expectation value through W = Z(T 2 — fxj). The 
new scale \xt can then be generated from dimensional transmutation, e.g., by replacing with 
a quark condensate as in Ref. [47]. For appropriate parameter choices, this extension leaves the 
O'Raifeartaigh dynamics largely intact. To connect the \xt scale to the A scale of supersymmetric 
QCD, we can introduce extra quarks Q' and Q' of SU(N C ) that make the theory conformal at 
high energies. The quarks obtain masses from W = TQ'Q', triggering the exit from a strongly- 
coupled fixed point. This provides an existence-proof model where all the scales are generated 
dynamically associated with single dimensional transmutation. The /i term can also be generated 
by introducing singlet "messengers" and coupling them, together with the doublet messengers, 
to the SSM Higgs fields as in Ref. [10]. 

Note that there is a Z2 parity in Eq. (149]) under which all the messenger fields are odd. Such 
a parity is often present in realistic constructions for the supersymmetry breaking sector. In 
order to avoid the problem of unwanted colored/charged relics, however, the messenger fields 
(either elementary or composite) must decay. If the decay occurs through dimension five oper- 
ators suppressed by M* and if the decay products contain an SSM state, then this may explain 
the discrepancy of the measured 7 Li abundance from the prediction of standard big-bang nucle- 
osynthesis, along the lines of Ref. [48J. In fact, the lowest dimension for the operators causing 
messenger decay can easily be five if, for example, the messenger fields are two-body bound states, 
such as meson states, of some strong dynamics. 

4 Astrophysical Signatures 

We have seen that dark matter in the present scenario naturally has the following features, which 
are not shared by the standard weakly interacting massive particle (WIMP) scenario: 
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• The mass of dark matter is of 0(10 TeV), which is significantly larger than the weak scale. 
The correct thermal abundance, however, is still obtained because of the relatively large 
annihilation cross section. 

• Dark matter can decay through dimension six operators, and thus with lifetime of O(10 26 sec). 
This can lead to observable cosmic ray signatures. 

• Dark matter can decay into light axion-like states (and to the gravitino if dark matter is a 
fermion), which in turn decays into standard model particles. The final states can naturally 
be only e + e~, 7r + 7r~7T°, or t + t~. 

In fact, these features are precisely what are needed to explain recent electron/positron cosmic 
ray data from PAMELA, FERMI, and H.E.S.S. In this section, we demonstrate that recent 
astrophysical data are indeed beautifully explained in the present setup, and provide a fit to the 
parameters m DM and r DM using the observed electron/positron fluxes. We also discuss the diffuse 
7-ray flux that could be seen in the near future by experiments such as FERMI. 

Let us begin with a summary of the observational situation. The PAMELA experiment has 
recently reported an unexpected rise in the positron fraction <3> e + / (<I> e + + $ e -) in the energy range 
between about 10 and 100 GeV [IB]. On the other hand, they did not see any deviation from the 
expected background in the antiproton data [IH] . The FERMI experiment is the first to be able to 
measure the combined electron and positron flux with good precision and control of uncertainties 
over the entire range from 20 GeV to 1 TeV. Their recent data also show an excess over standard 
diffuse cosmic ray backgrounds [21), with a broad structure extending up to the highest energiesS 
At higher energies, data from the H.E.S.S. experiment indicates a spectral break in the combined 
electron and positron flux at around 1 TeV, with the spectral index increasing from « 3.0 to 

« 4.1 [23]. 

While there remain experimental uncertainties as well as difficulty in calculating astrophysical 
background fluxes, taken together these results suggest a new source of primary electrons and 
positrons with a broad spectrum extending up to a few TeV. An exciting possibility is that 
these are signals of annihilation or decay of galactic dark matter J52J, [20] , although astrophysical 
interpretations, e.g. in terms of nearby pulsars [53], are also possible. When interpreted in terms 
of dark matter, the data has certain implications: 

• The structure around TeV in the combined e + and e~ flux is very broad, implying that 
electrons/positrons do not arise directly from dark matter annihilation or decay; rather, 
they arise through some cascading processes. 

14 The ATIC [50] and PPB-BETS [51] experiments have also reported an excess in the combined flux at around 
600 GeV, although with a peaked spectral shape that does not seem to fully agree with the FERMI data. With 
their significant statistical and experimental uncertainties, we do not include these results in our analysis. 
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• The spectral cutoff around TeV, together with the broad structure, implies that the mass of 
dark matter is larger than TeV. In particular, if cosmic rays arise from decay of dark matter, 
then the mass scale is more like O(10 TeV), especially if the cascade is sufficiently long. 

• The absence of signal in the antiproton data implies that dark matter annihilates or decays 
mainly into leptons, although precisely how much nucleonic final states should be suppressed 
is not completely clear because of uncertainties in proton/antiproton propagation models. 

These interpretations are consistent with detailed analyses performed after the recent FERMI 
data release [Ml 155] . 

The fact that a naive WIMP dark matter candidate, such as the neutralino lightest super- 
symmetric particle (LSP), does not satisfy the above criteria is quite suggestive. In the present 
scenario, dark matter mass is naturally of 0(10 TeV), which decays with lifetime of O(10 26 sec) 
producing the observed e signatures. The decay occurs through a light axion-like state (axion 
portal), so that the final states can selectively be leptons0 This also ensures that the final 
state leptons arise through cascades, making the spectrum consistent with the latest FERMI and 
H.E.S.S. data. It is also worth mentioning that decaying dark matter is much less constrained 
than annihilating dark matter [55], which has some tension with 7-ray and neutrino observa- 
tions [SHI E3- 

4.1 PAMELA, FERMI, and H.E.S.S. electron/positron data 

There are various possible decay chains in the present scenario, depending on which state is dark 
matter and the decay properties of the axion-like state a. Here we consider two representative 
classes: 1-step cascade ((f) — > aa; a — ► £ + £~) and 2-step cascade (0 — > 4>'4>'; 4>' — > aa; a — > 
where and represent dark matter and another state in the supersymmetry breaking sector, 
respectively, and i = e, //, r. We also consider a — > 7r + 7r~7T°. An example with a 1-step cascade 
was already seen in Section [3J dark matter is (the scalar component of) the lightest "meson," 
which decays into two axion-like states. The axion-like state then decays into standard model 
fields, as seen in Section 12. 5[ A 2-step cascade can easily arise if (the scalar component of) 
the lightest "baryon" is dark matter. After including dimension six baryon number violating 
operators, this state can decay into two meson states, each of which then decays into two axion- 
like states. Alternatively, dark matter may be a meson state which dominantly decays into other 
meson states, as seen in Section 13.31 

The processes just described are depicted in Figure [3J In the case where dark matter is a 
fermionic component, one side of some step of the cascade should be replaced by the gravitino, 

15 An R axion as the portal was mentioned in Ref. [18] and developed further in Rcf. [56 in the context of 
annihilating dark matter. For a related suggestion in dark matter annihilation, see [57] . 
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1 step 



2 step 



Figure 3: Cascade decays of dark matter through an axion-like state a. Here, <$' is an unstable 
state in the supersymmetry breaking sector, and £ (£ = e, /i, r) are standard model leptons. 

but this does not change the predicted signals, except that the dark matter lifetime must be 
rescaled by appropriate factors. Below, we present the fit of the predicted e ± fluxes to the 
PAMELA and FERMI data. We find that good agreements between the predictions and data 
are obtained for all the cases considered, except that the 1-step case may have some tension 
with the H.E.S.S. data. While we only present the results for a — > e + e~, and t + t~ here, 

we also performed the same analyses for a — > 7r + 7r~7r°. We find that the results in this case are 
very similar to the case of a — ► t + t~, since the final state spectrum arising from charged pion 
decay is similar to that arising from r decay. 

Our analysis here follows that of Ref. [SH], where galactic propagation of is treated by 
the standard diffusion-loss equation. The primary injection spectra are calculated as described 
there, assuming a large mass hierarchy in each cascade step. In addition, tau decays are simulated 
using Pythia 8.1 (60]. We assume an NFW halo profile [61] and use the MED propagation model 
parametrization given in [62] . While there is significant uncertainty in the halo profile within a few 
kpc of the galactic center, its effect on the predicted fluxes is small. We treat the astrophysical 
background as in Refs. [59l 52J: we take the parameterization of [63] and marginalize over the 
overall slope P and normalization A in the range —0.05 < P e +, e - < 0.05 and < A e + }€ - < 
oo. There remain substantial uncertainties in the background at the high energies explored by 
H.E.S.S., where primary electron fluxes depend strongly on individual sources within pa 1 kpc 
from the Earth. We therefore do not include the H.E.S.S. data in our fit procedure. 

We performed a y 2 analysis of signal plus background fluxes to the PAMELA and FERMI data. 
The PAMELA data at energies less than 10 GeV is strongly affected by solar modulation, and 
we exclude it from our analysis. The FERMI experiment released both statistical and systematic 
errors with their data. We conservatively combine these in quadrature, but note that this is 
likely an overestimation of the errors. The FERMI data is also subject to an overall systematic 
uncertainty in energy of "t^y, under which all energies are rescaled as E — > rE; we therefore also 
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Figure 4: Regions of best fit (at 68% C.L.) to the PAMELA and FERMf data for dark matter 
mass mDM and lifetime tdm, in the case of direct (solid), l-step (dashed), and 2-step (dotted) 
decays into e + e~, and r + r~. The best fit values of mDM and tdm are indicated by the 

crosses, and are displayed inset in units of TeV and 10 26 sec, respectively. Direct decays into e + e~ 
does not give a good fit. The case of n + n^n is similar to that of t + t~. 

marginalize over r in the range 0.9 < r < 1.05. The result for best fit regions of dark matter mass, 
771dm, and lifetime, tdm is shown in Figure IU for l-step and 2-step cascades to electrons, muons 
and taus. For comparison, we also show the fits for direct decays. With 26 degrees of freedom 
(7 PAMELA + 26 FERMI — 7 fitting parameters), we plot 68% CL contours, corresponding to 
X 2 = 28.8. The best-fit values of m D M and t D m are indicated in each case. 

We find that good fits are obtained in the region m D M ~ 0(1 - 100 TeV) and t D m ~ 
O(10 25 - 10 27 sec), depending on the decay chain. In Figs. [5] and [6] we show the predicted 
e ± fluxes compared to the PAMELA and FERMI data, using the best fit parameters obtained 
from the x 2 analysis. The H.E.S.S. high energy data [22] is also overlaid, with the energy rescaled 
in each case to best match the predicted flux, within the ±15% range of overall systematic un- 
certainty in H.E.S.S. energy. The agreement between the predictions and data is remarkable. We 
find some tension with the H.E.S.S. data in the case of l-step to e + e~, but in all other cases the 
predicted curves are consistent with the H.E.S.S. data despite the fact that we did not use it in 
our fit procedure. Note that the background fluxes are very uncertain above ~ 1 TeV, so that 
the precise comparison with the H.E.S.S. data may be misleading; for example, the background 
spectrum adopted here seems a bit too hard at the highest energies to be consistent with the 
H.E.S.S. data. 

To summarize, we find that the electron/positron fluxes observed by PAMELA and FERMI 
are very well explained by dark matter decay in our scenario. The mass of the axion-like state 
can take almost any value in the wide range considered in Section 12.3} except for a small window 
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Figure 5: The predicted fluxes compared to the PAMELA and FERMI data for 1-step cascade 
decays into e + e~, /i + yU~, and t + t~. In each case, the mass and lifetime of dark matter are chosen 
at the best fit point indicated in Figure HI with the background (dotted) and FERMI energy- 
normalization marginalized as described in the text. We overlay the H.E.S.S. data with energy 
rescaled in the range ±15% to best match the theory. Note that due to considerable uncertainty 
in the background fluxes at H.E.S.S. energies, direct comparison of predicted fluxes with the 
H.E.S.S. data may be misleading. 
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Figure 6: The same as Figure [5] for 2-step cascade decays into e + e , , and t + t 
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Figure 7: The high energy diffuse 7-ray spectrum away from the galactic plane, averaged over 
galactic latitudes above 10° assuming an NFW profile. Shown are the best fit parameters for 1- 
step to t + t~ (dashed), 1-step to 7r + 7r~7r° (dot-dashed), and 2-step to (dotted) dark matter 

decay modes. The 7-rays are due to ir° decay in the case of 1-step to t + t~ and 7r + 7r~7r°, and to 
FSR in the case of 2-step to fi + fi~. The r + r~ and 7r + 7r~7r° modes produce a bump in the flux 
clearly distinguishable from the background. 



at ~ 1.9 - 3.6 GeV where a decays hadronically. More specifically, we find that the regions 

(i) a^e+e- (2m e < m a < 2m M ) f a « O(100 TeV), 

(ii) a — > (2m M < m a <> m K - m^) f a « O(10 - 100 TeV), 

(ii) a -> //+//- (mjc - < m a ^ 800 MeV) / a w 0(1 - 100 TeV), (50) 

(iii) a -> tt+tt-tt (800 MeV < m a < 2m p ) / a w 0(1 - 100 TeV), 

(iv) a -> r+r- (2m T < m a < 2m fe ) / a w 0(1 - 100 TeV), 

can all explain the cosmic ray e ± data, without conflicting the bounds discussed in Section 12.51 



4.2 Diffuse gamma ray signals at FERMI 

An immediate consequence of the framework described here is that the decay of dark matter will 
provide a source of 7 rays throughout the dark matter halo, and extending up to energies around 
a few TeV. If the axion-like state decays into r + r~ or 7r + 7r~7r°, photons are produced directly by 
the decay of 7r°s, while for e + e~ and fi + fi~ modes there is a much smaller but significant source 
of 7 rays from final state radiation (FSR) and inverse Compton scattering (ICS). Although the 
greatest fluxes would originate from the galactic center, where number densities are highest, the 
best direction to look for them is away from the galactic plane, where the background is much 
smaller and the signal still large. The FERMI experiment will measure 7-ray fluxes over the entire 
sky at energies up to several hundred GeV, which has the potential to resolve spectral features 
caused by dark matter decays [6U [65] • 
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In Figure [7] we plot the contributions to the diffuse 7-ray fluxes for three illustrative decay 
modes: 1-step to r + r~, 1-step to 7r + 7r~7r°, and 2-step to yU + /i~. The first two cases are represen- 
tative of any cascade to taus or pions, while the fi + fi~ curve illustrates the much lower flux for 
a decay mode with only FSR photons. Assuming an NFW profile, we average over all galactic 
latitudes greater than 10° from the galactic plane, representative of a diffuse 7-ray measure- 
ment by FERMI. The astrophysical background is modeled by a power law flux, oc E~ 2 ' 7 , taken 
from [66J. For the fi + fi~ mode the flux shown is solely due to FSR, and we assume an a mass 
of 600 MeV. We have not shown the contribution from ICS, which, like FSR, is subdominant to 
directly produce photons for t + t~ and 7r + 7r~7T° modes. At its peak, the ICS flux- component is 
expected to be comparable to the background, with a spectral shape somewhat different to that 
from FSR [551 [65], an d could be significant for e + e~ and fi + n~ modes. 

We find that signals from dark matter in our scenario can be seen in the diffuse 7-ray data. 
While the background is very uncertain, it is expected to be smooth compared to the strong, 
peaked fluxes seen for the r + r~ and 7r + 7r~7r° decay modes, which should result in a clearly visible 
bump in the spectrum. The weaker signals from e + e~ and fi + fi~ modes may also be seen. Since 
astrophysical sources are not expected to produce 7-ray fluxes with such prominent spectral 
features, measurements of diffuse 7-rays may serve to distinguish dark matter as the source of the 
PAMELA and FERMI excesses. Additionally, the shape and the strength of the 7-ray spectrum 
would convey information about the mass and decay channels of the dark matter. While the 
absence of an excess in upcoming FERMI data would not exclude our scenario, especially for the 
e + e~ and cases, a positive result would be a striking signature. 

5 Collider Signatures 

If the PAMELA and FERMI data is indeed indicative of dark matter with an O(10 TeV) mass, 
then direct production of dark matter is not possible at current collider experiments. However, 
the light axion-like states, which play a crucial role in determining the cosmic ray spectra, are 
kinematically accessible. We briefly describe some of the collider signatures for the axion-like 
state a. 

The discovery potential for the axion-like state depends strongly on its mass and decay con- 
stant. Since we are considering decay constants in the range f a ~ 0(1 - 100 TeV), the coupling 
of a to the standard model is small, so discovery relies on the axion-like state having a clean 
decay mode. For 2m e < m a < 2m M , the decay length is generically greater than a kilometer, 
so any axion-like state produced in a collision will decay outside of the detector. In mass range 
2m fl < m a ^ 800 MeV, a decays into /U + /i~ with a possible displaced vertex, offering a promising 
discovery channel. For 2m T < m a < 2mf,, a decays promptly into taus, but neither leptonic taus 
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nor hadronic taus are particularly clean channels. Similarly, there are large hadronic backgrounds 
to a — > 7r + 7T~7r°. Interestingly, while the r + r~ and 7r + 7r~7r° channels are challenging in the col- 
lider context, they are precisely the ones that give the most striking diffuse 7-ray signal as seen in 
Section H~2l In contrast, the fi + n~ final state yields less dramatic diffuse 7-ray flux, but relatively 
clean collider signatures. 

Direct production of a's at the LHC was considered in Ref. [21], where the a is produced 
in association with a hard (pt > 100 GeV) jet via the gluon-gluon-a coupling. Since the a 
is boosted from the production, the resulting muon tracks from a decay have a mrad opening 
angle. If the decay constant /„ ^ (3 - 10) TeV, then the axion-like state lives sufficiently long 
that the decay a — > \T happens with an O(cm) displaced vertex. As long as the decay 
constant /„ ^ (15 - 30) TeV, then the production rate is sufficiently large to compete with the 
B° — > D /j^v — > ^Tvv background. Therefore, direct production of a is a promising possibility 
for f a ~ 10 TeV. 

An alternative production mechanism for a's is via the Higgs boson. For the R and PQ sym- 
metries considered in Section 12.31 the axion-like state has couplings to the Higgs fields, allowing 
the Higgs boson to decay into two axion-like states. For simplicity, we will focus on the decoupling 
limit as in Eq. (122]) . The Higgs decay width into two a's is then given by 

r( " -* aa) = ikrjt' (51) 

where c\ = (1/4) sin 2 2/?, and rrih is the Higgs boson mass. For a light Higgs boson, the dominant 
decay mode is into bb with (taking m^jm^ — > 0) 

T(h^bb) = ^-m h (^-) 2 . (52) 

Therefore, the branching ratio into two axion-like states is 

s c?m?t» 4 , / m h \ 2 /i0TeV\ 4 

where we have used the large tan j3 approximation, sin2/3 « 2/ tan (3, in the last equation. We 
see that the branching fraction for h — >• aa depends strongly on the decay constant as f~ 4 , and 
a sizable branching fraction can be obtained as f a approaches 1 TeV. Since the decay constant 
is naturally expected to be somewhat (« Air) smaller than the dynamical scale A as we saw 
in Eq. (j!4j) . this makes observation of the h — > aa — > 4/i signal at the LHC an interesting 
possibility [IB]. The recent analysis of Ref. [SZ] focused on the h —>■ aa — > 2/i2r channel in 
the case of m a ~ 5 GeV. They find that generic event selection has an efficiency around 10%. 
Assuming a similar efficiency for h — > aa — > 4/i events, and taking the Higgs production at the 
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LHC of ~ 50 pb, at least 10 events could be seen with 300 fb" 1 of data for f a tan (3 £ 10 TeV. In 
this way, direct a production (for larger values of f a ) and a production through the Higgs boson 
(for smaller values of f a ) are complementary. For larger m a , the h — > aa — > 6tt and h — > aa — > 4r 
channels may be visible if the background can be controlled. Note that the recent D0 analysis 
in Ref. [68] already gives the constraint / a tan/3 2 TeV for a — > 

There has also been recent interest in looking for light bosons in low-energy high-luminosity 
lepton colliders [69j EO] as well as in fixed-target experiments [70]. From Eq. (fl9"]) ; the coupling 
of the axion-like state to electrons is proportional to m e / f a , which is smaller than 10~ 6 for 
f a > 1 TeV, making the process e + e~ — > 7a beyond the reach of current lepton colliders. The 
feasibility of a fixed-target experiment to discover a depends on the a lifetime. If a decays 
promptly, then one must contend with a huge standard model background of prompt charged 
particle production, though Ref. [70] suggests that a coupling as small as 10~ 6 might be testable 
if a decays into If a decays with a displaced vertex, then it could be discovered in traditional 

electron or proton beam dump experiments. Given the bounds from Figure [2] and the lifetimes 
in Eqs. < H7\ I2B"1 [251 , the possible values of cr spans a huge range from tens of kilometers to less 
than a nanometer. If a — > /i + yU~, then cr is plausibly in the millimeter to centimeter range, and 
could likely be discovered in an upgraded version of the experiment from Ref. [38] . 

6 Discussion and Conclusions 

The origin of dark matter is one of the greatest mysteries in particle physics and cosmology. From 
the theoretical point of view, attentions have been focused on the WIMP paradigm: dark matter 
has mass ~ O(100 GeV - 1 TeV) and couplings of weak interaction strength g m 0(1), 

leading to an annihilation cross section that gives the correct thermal relic abundance, ^dm — 0.2. 
However, since the annihilation cross section depends on the combination g 2 / m DM, it should be 
equally convincing to consider the case where tit-dm is heavier than the weak scale as long as the 
coupling g is larger. Such a situation arises naturally if dark matter is a composite state of some 
strong interaction, since the typical coupling g is expected to be larger than order unity and ttidm 
can be O(10 TeV). 

Suppose that dark matter indeed arises from some strongly interacting sector. Then its stabil- 
ity may be the result of compositeness, and not of some exact symmetry imposed on the theory. 
This is precisely analogous to the case of the proton in the standard model embedded in some 
unified theory. If the proton (and pions) were elementary, it would immediately decay into e + 
and 7T° through a Lagrangian term £ ~ pen . Since it is composite, however, the leading operator 
causing proton decay is already dimension six, £ ~ qqql/M^, and the resulting lifetime is of order 
10 36 years for M* of order the unification scale. In the case of composite dark matter, dimension 
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six operators yield a lifetime of order 10 25 sec for /tidm ~ 10 TeV and M* « 10 GeV. Decay of 
galactic dark matter could then have currently observable consequences. 

The story just described implies that there is new strongly coupled physics beyond the weak 
scale, at « O(10 - 100 TeV). Interestingly, we already know an attractive framework where 
such a situation occurs — weak scale supersymmetry with low energy dynamical supersymmetry 
breaking. Since the superparticle masses in this framework arise at loop level, the scale of the 
strong sector is naturally larger than the weak scale by a one-loop factor « 16ir 2 . This picture is 
very much consistent with what is implied by the LEP precision electroweak data, namely that 
physics at the weak scale itself is weakly coupled. Yet such a setup can still explain the large 
hierarchy between the weak and Planck scales in, arguably, the simplest manner via dimensional 
transmutation. 

It is interesting that this picture of dark matter arises precisely in the scenario where conven- 
tional LSP dark matter is lost — the LSP is now the very light gravitino. In fact, a gravitino with 
mass m 3 / 2 ^ O(10 eV), implied by A « O(10 - 100 TeV), avoids various cosmological difficulties 
faced by other supersymmetry breaking scenarios. This allows us to consider the "standard" cos- 
mological paradigm, based on inflation at a very early epoch with subsequent baryogenesis at high 
energies, consistently with supersymmetry. The standard virtues of weak scale supersymmetry, 
such as the stability of the weak scale and gauge coupling unification, are all preserved. 

While composite dark matter in low-scale supersymmetry breaking already offers a consistent 
picture, this may not be the end of the story. Since dark matter is a part of a strongly interacting 
sector, it is possible that it feels additional dynamical effects. In particular, it is quite plausible 
that the sector spontaneously breaks an accidental global symmetry, leading to a light pseudo 
Nambu-Goldstone boson — again as in QCD. This raises the possibility that dark matter decays 
mainly into these light states (possibly through some other state), which then decay into standard 
model particles. In supersymmetric theories there is a natural candidate for such symmetry, 
an R symmetry, whose existence is suggested by a certain genericity argument associated with 
supersymmetry breaking. The mass of these light states can easily be in the MeV to 10 GeV 
range if explicit breaking arises from dimension five operators. Except in the mass range m a ~ 
1.9 - 3.6 GeV, the decay products of the light states are naturally "leptonic" (specifically, e + e~, 
/U+yU - , 7r + 7T~7r°, or r + r~), with little nucleonic activity. 

An illustrative model of this class was given in Section [3j and the required symmetry structure 
for the supersymmetry breaking sector was summarized in Table [TJ Remarkably, the properties 
we have just obtained (with a little adjustment of parameters), are precisely what is needed to 
explain the recent cosmic ray data through dark matter physics: 

• Dark matter mass is of O(10 TeV), 

• Dark matter lifetime is of O(10 26 sec), 
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• Dark matter decays through (long) cascades, 

• Dark matter decay final states are leptonic. 

As we saw in Section HI various mysterious features in the data — an unexpected rise of the 
positron fraction between w 10 and 100 GeV in the PAMELA data, nonobservation of any 
anomaly in the PAMELA antiproton data, and a broad excess of the combined e ± flux in the 
FERMI data — are all beautifully explained by the properties of dark matter discussed in this 
paper. The resulting spectra are also consistent with the recent H.E.S.S. result. The success 
is quite remarkable, especially in view of the fact that the data is difficult to explain in terms of 
conventional WIMP annihilation — the mass scale suggested does not seem natural, the observed 
rate requires a large boost factor, and typical WIMP annihilation produces more antiprotons than 
indicated by PAMELA. 

While dark matter with mass of order 10 TeV will not be produced at the LHC, the present 
scenario still has potential collider signatures. Since the light states generically have interactions 
with standard model gauge bosons, they may be produced at the LHC, and the leptonic decay of 
the light state could lead to visible signatures, especially if the dominant decay is into muons. The 
Higgs boson may also decay into the light states, producing a clean four lepton final state. Future 
astrophysical observations could also probe this scenario. For example, measurements of diffuse 
7-ray could discriminate e//x final states from t/tt final states. Moreover, the present framework 
may explain the discrepancy between the measured 7 Li abundance and the standard big-bang 
nucleosynthesis prediction by dimension five decay of some of the states in the supersymmetry 
breaking sector, e.g. (composite) messenger fields. Analysis of all these experimental data could 
provide important information about the structure of the supersymmetry breaking sector. 

We finally mention that while in this paper we focused on the case where dark matter arises 
from a strongly interacting supersymmetry breaking sector and decays through light states, some 
of our results apply in wider contexts. For example, quasi-stable dark matter in the supersym- 
metry breaking sector may decay directly into the SSM sector particles through dimension six 
operators. For example, quasi-stable mesons can decay into SSM (s)leptons Lssm through the 
Kahler potential operators K ~ Q[Q^ L, ssu Lssm/ M% + Q^Qj-^ssm-^ssm/M? . In these cases, the 
final states of the decay are determined by a combination of gravitational scale physics and TeV- 
scale superparticle spectra, not through cascades associated with light states. However, all other 
points regarding compositeness, long lifetime, relatively large mass, and thermal abundance still 
persist. The dynamics of dark matter discussed here can also be applied in non-supersymmetric 
theories — all we need is some strong dynamics at « O(10 - 100 TeV) satisfying the properties of 
Table [TJ Examples of such theories may include ones in which the Higgs fields arise as pseudo 
Nambu-Goldstone bosons of some strong dynamics [71] . 

As the LHC starts running this year, we anticipate great discoveries. A possible connection 
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between weak scale physics and dark matter is one of the major themes to be explored at the high 
energy frontier. The standard expectation is to study the properties of a dark matter particle 
by producing it at the LHC. Although the scenario presented in this paper does not allow this, 
important physics associated with the dark sector may still be probed. With many new particle 
physics and astrophysics observations on the horizon, the next decade will certainly be exciting for 
potentially understanding the origin and properties of dark matter. While nature may not show 
us the "standard" WIMP dark matter story, it may still reveal a beautiful connection between 
dark matter and the weak scale through hidden sector dynamics. 
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